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Inspecting absorption in the spectra of extra- galactic gamma-ray 
sources for insight into Lorentz invariance violation 

Uri Jacoqj and Tsvi Piran 

Racah Institute of Physics, The Hebrew University, Jerusalem, Israel 

We examine what the absorbed spectra of extra-galactic TeV gamma-ray sources, 
such as blazars, would look like in the presence of Lorentz invariance violation (LIV). 
Pair-production with the extra-galactic background light modifies the observed spec- 
tra of such sources, and we show that a violation of Lorentz invariance would gener- 
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■^j- . ically have a dramatic effect on this absorption feature. Inspecting this effect, an 



experimental task likely practical in the near future, can provide unique insight on 
the possibility of LIV. 



INTRODUCTION 



(N 
> 

Quantum-gravity theories predict in general the breakdown of familiar physics when 
approaching the Planck energy scale (E p i = ^Jhc b /G ~ 1.2 x 10 28 eV). A violation of 
Lorentz invariance arises naturally in various such theories (for a review see e.g. 

HQ). 

Though the possible effects of LIV are not currently measurable on Earth, there is vast 

O ■ 

interest in exploring LIV, as it could direct us to the fundamental physical theory. We 
would like to obtain experimental evidence for a LIV scenario or, to the contrary, evidence 
strengthening the validity of Lorentz invariance and constraining the possibility of LIV. 

We consider here a general approach that involves the breakdown of Lorentz symmetry 
at some very high energy scale, £E p i. At lower energy scales small deviations from Lorentz 
invariance may be felt. This kind of LIV approach is inspired by the theoretical work 
on quantum-gravity, but the phenomenology may of course be considered independently of 
quantum-gravity. The basic phenomenology of LIV can be described by a modified dispersion 
relation. When examining particles with E -C £E p i, we may regard only the leading order 
correction to Lorentz invariance and generically (for ultra-relativistic particles) write the 
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dispersion relation as: 

E * -j?c 2 -m 2 c 4 ~±E 2 (-^- \ , (1) 



where the dimensionless parameter, £ n , and the order of the leading correction, n, are model 
parameters (which can be universal or be given separately for each particle type, depending 
on the framework), c, the conventional speed of light constant, remains only the speed of 
low-energy massless particles, and we put hereafter c = 1. We notice, assuming that the 
standard relation v = dE/dp holds, that the + and — signs account for energy- dependent 
super luminal and infraluminal motions. 

This dispersion relation points us to the most direct method of testing LIV - inspecting 
particles' speeds. As the effects of LIV may become substantial for high-energy particles, 
there has been much work recently on testing LIV by means of time-of-flight measurements 
of astronomical particles. Time-of-flight analyses provide the most generic method of testing 
LIV and have been performed for different cosmological gamma-ray sources. No signatures 
of LIV have been detected, yet these observations were at most sensitive to effects of LIV 
with a scale two orders of magnitude below the Planck scale [siO,^, 6, 3, 2, 9|. Time-of-flight 
measurements of cosmological neutrinos are expected to be more sensitive 10]. 

Another observational option for examining LIV is to inspect the thresholds of high- 
energy particle reactions, as these might be shifted as a result of a LIV alteration of the 
transformations between reference-frames. This test is not strictly kinematic and so less 
generic. Nevertheless, the thresholds can be derived directly from the deformed dispersion 
relation (with the sole assumption being that it represents the only modification of standard 
physics). Regardless of the specific theoretical details and the particle dynamics within 
a given framework, one can expect that theories where the Lorentz symmetry is broken 
at extreme energies would be manifested in the low-energy regime by an approximated 
dispersion relation of the form ([I]). Thus, this method for gaining insight on LIV, dealing 
with particle interactions, is not model-specific and is also an efficient approach. 

The cutoff in the spectrum of ultra high energy cosmic rays, the so-called GZK effect 
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121 ]. which in the past failed to be observed as expected, created an apparent paradox 
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LIV 



and originally motivated the consideration of LIV as a possible solution |13| 
could have pushed upward the threshold energy of the reaction between the cosmic rays 
and the cosmic microwave background (CMB) photons. The GZK cutoff has more recently 
been detected by the HiRes observatory [l6| and then by Auger with high significance 
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171 ]. Turning the previous argument around, the agreement of the GZK threshold with 



the classical expectation can be used [18(] to place a very severe bound on LIV well above 
the Planck scale. However, some theoretical frameworks predict a priori that LIV does not 
affect cosmic rays (as any charged particles), making the above bound less generally valid 



ial2Q|,l21] 



Considering that in some LIV frameworks only photons deviate from classical behavior, 



phenomena invo 
suggested 
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just photons are the most suitable LIV tests. Specifically, it has been 



24j that LIV could affect the cosmological interaction between gamma- 



rays and infra-red background photons. We proceed to demonstrate here the feasibility of 
testing LIV by investigating absorption in the TeV spectra of gamma-ray sources. We show 
that this can supply unique evidence for LIV (or rule out its existence) at the Planck scale 
level. 

The paper is organized as follows. In Sec. [TT] we review the effect of LIV on the pair- 
production interaction threshold. Section [Till describes the extra-galactic background field 
and its impact on the optical depth and observation of gamma-ray sources. In Sec. IIVI we 
explain and demonstrate how LIV would affect the absorption features of cosmological TeV 
gamma sources. We discuss our results and the practicality of gaining valuable insight on 
LIV in Sec. ED 



II. DEFORMED INTERACTION THRESHOLDS 



Numerous particle processes, such as the two mentioned above, only take place beyond 
a threshold, required for the conservation of energy-momentum and classically determined 
by Lorentz-transforming to the center-of-mass reference frame. The photopion production, 
responsible for the GZK effect, does not occur if the energy in the center-of-mass frame is not 
sufficient for the creation of a pion. The pair creation interaction, in which we are interested, 
can happen if in the center-of-mass frame the photons have a total energy larger than the 
rest mass of the electron-positron pair. Within LIV the transformation to the center-of-mass 
frame is deformed and the threshold energy in the laboratory frame can be shifted. Since in 
the context of a generic LIV approach we do not have global laws for transformation between 
reference frames, we must work fully in the laboratory frame. The interaction threshold can 
still be determined by equating the energy and momentum of the particles before and after 
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the reaction, using the particles' dispersion relations and requiring that the solution to the 



equations is real at the threshold or above and imaginary below it (see 15|). The single 
assumption made here is that standard energy-momentum conservation is maintained in 
oarticle reactions (note that this assumption does not hold in Doubly Relativity models 



|2* 



The classical prediction of Special Relativity is that the threshold energy of a soft photon 
for pair-production by head-on collision with a gamma-ray of energy E 1 is m^/E^, where 
m e is the electron rest mass. Within our LIV phenomenology we use ([1]) and arrive, as 
explained, at the modified pair-production threshold: 

e thr = mj/E, + I {£^j E 7 , (2) 

if the deformed dispersion relation is only applied to photons, or 

e thr = ml/K, + ~(1 - 2"») (£^j " E„ (3) 

if LIV affects all particles equally. We want our results to be relevant for theories where the 
Lorentz- violating physics applies only to photons (and not electrons) as well as for theories 
where the LIV correction is universal (natural when the correction is due to the small-scale 
structure of spacetime and the equivalence principle holds). However, we find, as far as 
the interaction of gamma-rays with background photons is concerned, a similar LIV effect. 
The application of LIV to electrons accounts for the reduction of the LIV term in ([3]) by 
a maximal factor of 2. This is equivalent to modifying £ by such a factor, and we show 
subsequently that our results depend very weakly on the choice of the symmetry breaking 
scale. We, therefore, use for the rest of our study the threshold energy expression of 
as a more conservative estimate of the LIV effect, and our findings will be generically valid. 
The threshold expression of ([3]) includes the leading order correction in the small factors 
E Z and ^ and accounts for a LIV scenario with subluminal motion [' — ' sign in Eq. (fT])l. 
Such a LIV scenario, which produces an upward shift of interaction threshold energies, is 
more probable theoretically and observationally, and we proceed to study it. 1 



1 A LIV scenario with superluminal motion would have the correction term in ([3]) negative. The interaction 
threshold energy would be lowered, approaching zero above E*. This would imply that gamma-rays of 
energy ~ E* can interact with the bright CMB field, resulting in a cutoff of all extra-galactic gamma-ray 
sources around E* . This LIV scenario therefore has a very clear signature, and it is unlikely with current 
observations. 
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Inspecting ([3]) we find that in the context of LIV there is a critical energy, which we 
denote E*, for which e t h r becomes minimal. This new phenomenon at E* is a key ingredient 
in our approach. The critical energy is straightforwardly derived by setting e' thr (E 7 ) = 0. 
The general expression is: 



It is apparent that the dependence on the symmetry breaking scale is weak, especially for 
n — 1,2. With n = 1 varying £i by two orders of magnitude results only in changing E* 
from 11 TeV for £i = 0.1 to 50 TeV for £i = 10. We happen to be lucky that this is a 
convenient energy range and for £i = 1 (the symmetry breaking scale being the Planck 
energy) the critical E* is at the examinable energy of 23 TeV. For n = 2 E* ~ 10 5 TeV 
with the symmetry breaking around the Planck scale. Scenarios with n ^ 2 are, therefore, 
clearly unapproachable. We carry on demonstrating the applicability of our method for the 
investigation of the theoretically motivated n = 1, £ ~ 1 LIV scenario. 

We present in Fig. [1] the shift of pair-production threshold energies caused by LIV. It is 
apparent that at low gamma-ray energies the effect of LIV is negligible, while it becomes 
qualitatively significant when approaching E*. For E* the soft photon threshold energy 
reaches a global minimum and then increases with increasing gamma-ray energy (as opposed 
to the systematic decrease in the Lorentz invariant case). We notice that some soft photons 
(at energies below ~ 0.02 eV with the LIV parameters here) never interact with any hard 
photons, as they are below the pair-production threshold for every gamma-ray energy. For 
every soft photon that does interact with gamma-rays above a lower threshold, there is 
also an upper gamma energy threshold above which the interaction does not take place. We 
explore the implications of these peculiar high-energy interaction properties in the following. 

III. THE EXTRA-GALACTIC BACKGROUND RADIATION AND 
COSMOLOGICAL GAMMA-RAY SOURCES 

The extra-galactic medium contains uniform background radiation fields in several spec- 
tral ranges. The part of the spectrum that consists of the infra-red, optical and ultra-violet 
backgrounds is commonly referred to as the extra-galactic background light (EBL). The 
EBL, which shines at energies just above the familiar CMB, is thought to consist of the 
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FIG. 1: The soft photon pair-production threshold energy as a function of the interacting hard 
photon energy - according to the classical calculation of Special Relativity (dashed black line) and 
with n = 1, £ = 1 LIV (blue curve). 



accumulated redshifted thermal emission of stars and dust throughout the history of the 
universe. Yet, as opposed to the well-resolved CMB, this background is difficult to de- 
termine, as direct measurements are problematic, owing to the bright galactic and solar 
system foregrounds present. Certain estimates and limits on the EBL spectrum exist. In 
the infra-red range, of highest relevance to our study, they arise mostly from COBE (Cosmic 
Background Explorer) and IRTS (Infrared Telescope in Space) observations and from galaxy 



counts. Still, there are significant uncertainties regarding the EB 



26 



spectrum, and theoretical 
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29). 



estimates are very cosmological-evolution model-dependent 

Gamma-rays from cosmological sources are attenuated by pair-production interaction 
with the EBL, so a source appears on Earth with a different spectrum than its emitted 
one. The observed gamma-ray spectrum from a source at redshift z s is given by I b s [E^] = 
Iint[E<y(l + z s )\e~ T ( El,Za \ where Im is the intrinsic source spectrum and r is the respective 
optical depth. The calculation of the optical depth uses the pair-production cross-section 



30t ] and is determined by the spectral energy distribution of the EBL and the length of the 
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inter-galactic path over which the gamma-rays propagate: 



dl(z) / de (T 77 [E 7 (l + z),e(l + z))n{e)(l + zf 



(5) 



where dl is the proper line element (we use standard cosmology parameters), <r 77 is the 
cross-section for 77 —> e~e + in an isotropic field of background photons and n(e) is the EBL 
photon number distribution at present. We employ a simplistic model for the cosmo logical 
evolution, where n^[e(l + z)] = n(e)(l + z). This is negligible in our study which is concerned 
with low redshift sources. 

Over recent years several extra-galactic sources (mostly blazars) have been observed by 
imaging air Cherenkov telescopes (IACTs) in the TeV gamma-ray range. The spectra of 
some observed blazars seem to be harder than anticipated, considering the expected EBL 
absorption. However, this expectation is ambiguous because of the uncertainty in the EBL 
spectrum. Several studies have, thus, attempted to extract a more reasonable EBL spectrum 
out of the observed blazars spectra. This is also somewhat problematic, since our knowledge 



of the intrinsic sources spectra is limited 
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34( . Other studies have 



that the harder than anticipated spectra could be due to LIV effects [15l . l24j . Given the 



urther suggested 



uncertainties it seems unreasonable to justify an exotic explanation such as LIV effects in the 
context of such studies. We use the exceptional structure of the Lorentz- violating interaction 
threshold to suggest in this study a simple analysis of observational data that can give a 
clear sign of LIV or conclusively rule it out regardless of EBL specifics. 



IV. ABSORBED GAMMA-RAY SOURCES IN THE PRESENCE OF LIV 

The pair-production reaction's cross-section is concentrated in a narrow energy range 
above the threshold (peaking at ~ 4 times the threshold energy). As classically the soft 
photon interaction threshold is a monotonous function of the gamma-ray energy, given that 
lower energy background photons are more abundant than higher energy ones [the spectrum 
is steeper than n(e) = e _1 in the relevant range], the optical depth becomes monotonously 
larger for higher energy gamma-rays. For every cosmological source distance there is a 
gamma-ray cutoff energy at which the optical depth becomes unity, and around it the 
observed spectrum should show exponential suppression. LIV would change this behavior 
drastically. In the presence of LIV there is a certain gamma-ray energy, E* (see Sec. HT|) . 
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above which higher energy gamma-rays require higher energy EBL photons for the pair- 
production interaction. For any cosmological source the optical depth would thus be maximal 
at E*, and at higher energies it would rapidly decrease. 

In the classical calculation of the optical depth (see 3(3]) the pair-production cross-section 
is determined by the soft photon's energy relative to the threshold for interaction with a hard 
photon of E 7 , i.e. a^[E^,e] = cr 7T [e/e t ^ r (£' 7 )]. To proceed within a generic LIV approach 
we assume that the modified cross-section, a 77 £jy, has the same dependence on e/e t hr- The 
cross-section is merely shifted as we now evaluate it by regarding the soft photon's energy 
relative to the LIV shifted pair-production threshold - we use the e t hr of (EJ) to evaluate 
e/e thr (E 7 ). When e t hr,Liv{E) = €thr, classical 

(E ), we have ct 7% liv [E, e] 

&fy t classical [E ,e], 

and we can put in (jSJ) cr 17t Liv instead of a^, classical ■ This is natural within any relativistic 
framework, but since the LIV framework might have a preferred inertial frame, the implied 
assumption is that modifications to the shape of the pair-production cross-section would 
constitute lower order corrections. With this basic assumption we find at any energy above 
E* an optical depth identical to the optical depth at some energy below E* for which the 
pair-production threshold is the same. This "symmetric" behavior is the key to our proposed 
LIV test. 

The effect of the modified pair-production threshold on the optical depth of gamma- 
rays can be seen in Fig. [2j We used for simplicity an EBL spectrum of n(e) = 10~ 3 e~ 2 
[eV -1 cm~ 3 ] (approximately the average energy dependence of the EBL at the relevant infra- 
red energies) , as the LIV features examined are not sensitive to this choice. We notice that 
for sources at any cosmological distance the optical depth is maximal at E* = 23 TeV. For 
some sources the optical depth reaches unity at a lower energy, increases until 23 TeV and 
then decreases back to unity and below at higher energies. More nearby sources have optical 
depths that never reach unity, as they have r(23TeV) < 1. While the detailed shape of the 
optical depth at gamma-ray energies up to E* depends on the detailed densities of the EBL 
spectrum, the optical depth above E* does not depend on additional background spectrum 
data. At any energy above E* the optical depth is the product of interaction with the same 
part of the EBL spectrum that produces the optical depth at some energy below E* (see 
Fig. [1]). The correspondence between values on both sides of E* is strictly dependent on the 
LIV parameters. This means that whatever the details of the EBL and the optical depth 
below E* are, with the LIV model we can foretell what the optical depth above E* will look 
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FIG. 2: The optical depth for extra-galactic sources at sample redshifts as a function of gamma-ray 
energy. We use a simplified EBL spectrum (see text) and present the results with n = 1, £ = 1 
LIV (blue curves) and without LIV (dashed black lines). 



like. For every source redshift from which a cutoff of the gamma-ray emission by the EBL 
absorption is observed, we can foretell the energy of re-emergence, where the optical depth 
decreases below unity again. 

Figure [3] shows how LIV would affect a blazar's observed spectrum. We demonstrate this 
with a more realistic realization of the EBL spectrum, vl u = 10 nW/m 2 /sr below 1.5/im, 
going as A -1 up to 10/im and as A 1 above (the EBL wavelengths relevant for attenuation 
of TeV gamma-rays are ~ 1-100/im). Of course this is yet one realization out of many 
possible (see 28|] ) . Our intention is simply to display the manifestation of the general 
effect. The figure also depicts the observed blazar's spectrum anticipated when considering 
the simplified EBL spectrum presented earlier. We notice that the principal LIV effects 
are the same. The sample blazar's intrinsic spectrum is chosen to be a power-law with 
a spectral index (dN/dE oc E~ a ) of a = 2. This is a likely approximate spectral index 



for blazar emission in the TeV range 



351 ] . We may expect to be able to estimate the 



intrinsic spectrum of a TeV gamma-ray source from lower energy measurements, where the 
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FIG. 3: A sample blazar's spectrum on Earth in the presence of LIV at the Planck scale. We 
assume a source at z = 0.03 with a simple intrinsic spectrum (see text) marked by the thick 
red line. The observed blazar spectrum is presented with the simplified EBL spectrum of Fig. 
[2] (dashed curves) as well as with a more realistic EBL spectrum (solid curves, see text). Blue 
curves describe the extra-galactic absorption within LIV, whereas black curves describe classically 
expected absorption. 

effect of the extra-galactic attenuation is negligible. Of course the source's intrinsic spectral 
shape may vary and there may be an intrinsic break in the spectrum, but the accumulation 
of observations of different TeV gamma-ray sources will help us separate intrinsic source 
features from those induced by the EBL interaction (see discussion later). 

We are also faced with the consideration of whether enough photons will be detected 
by the gamma-ray observatories at energies above E*, so that a source's spectral shape at 
these energies can be determined. Fortunately, the LIV scenario allows us in principle to 
observe nearby blazars at much higher multi-TeV energies than in the classical the 
optical depth decreases instead of rapidly increasing. We demonstrate this in Fig. H]with the 
expected photon fluxes from sources with different intrinsic spectra. Mkn-501, a powerful 
blazar at a redshift z = 0.034, has been observed up to ~ 20 TeV and was found to have 



11 



dN/dE (a.u.) 



0.01 



10- 4 



10" 



10" 



lxlO 12 2xl0 12 



5xl0 12 lxlO 13 2xl0 1; 



5x10° lxlO 1 



By (eV) 




FIG. 4: The expected photon fluxes from different sample blazars in the presence of LIV at the 
Planck scale. We examine three sources at a redshift z = 0.03 with intrinsic power-law spectra, 
dN/dEi n t oc E~ a . The blue solid curves describe the observed spectra of the sources with a = 
2,2.25,2.5 (from top to bottom), considering LIV absorption with the same EBL realization as 
the solid curves in Fig. [3j The dashed black curve describes the observed spectrum of the a = 2 
blazar with no LIV. The dashed red line corresponds to the intrinsic spectrum dN/dE oc E~ 2 . The 
normalized time-averaged fluxes of Mkn-501, measured during the 1997 outburst [36], are presented 
as black points with statistical errors. At higher energies we have 2a upper flux limits (systematic 
errors are not included). 



a spectrum that can be interpreted as a power-law of a ~ 2 with an exponential cutoff 
due to EBL extinction [3^]. As the n = 1, £ = 1 LIV model does not have a significant 
effect on the optical depth below 20 TeV, this observed spectrum can also be easily fit to 
absorption with LIV (where we may also use the freedom of the EBL realization). Above 
~ 20 TeV, however, LIV would qualitatively change the observed spectrum. The HEGRA 
observation of Mkn-501 during the extraordinarily active period in 1997 produced upper 
flux limits up to ~ 40 TeV (see Fig. H]). These upper limits are not stringent enough to rule 
out the LIV scenario - they seem to exclude an intrinsic power-law of a < 2 for Mkn-501, 
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but the data may still be reasonably fit to LIV with a ~ 2.25 (taking into account the 
systematic error arising from the 15% uncertainty on the absolute energy scale improves 
this fit). However, as apparent from Fig. HI we are on the verge of obtaining significant 
insight on LIV with symmetry breaking at the Planck scale. An observation of Mkn-501 
at slightly higher energies would clearly distinguish between a LIV and a no-LIV scenario. 
Stronger upper limits at the high energies would tend to rule out LIV. Such observations 
could be possible with improved observatories or with another high-activity period of the 
source. 

We see in Fig. 0] that while in the classical scenario the extra-galactic extinction becomes 
so strong that we cannot observe gamma rays much above 20 TeV even with a flatter intrinsic 
spectrum, within the LIV scenario it is likely that we will be able to detect photons at higher 
energies regardless of the exact intrinsic spectral index. As an example, if such a nearby 
blazar's intrinsic TeV spectrum is dN/dE oc E~ 2 , then, given LIV, the number of photons 
expected to arrive from it at the range of 20 — 40 TeV is only a factor ~ 3 smaller than the 
number at 10 — 20 TeV. Furthermore, at 40 — 80 TeV the number of expected photons is 
a factor ~ 2 smaller than at 10 — 20 TeV. Hence, if the LIV model holds, it is a realistic 
expectation to detect its signature with future observations. Figure H] also shows that this 
signature is quite independent of the intrinsic power-law shape. All in all, although we do 
not know the specifics of the blazars' intrinsic spectra or the EBL spectrum, we see that 
LIV would have such a dramatic and unique effect on the observed gamma-ray spectra that 
we can expect to notice it and single it out from other effects. 

Having demonstrated the detectability of the LIV effect for the scenario of £i = 1, we 
should note that if the symmetry breaking scale is much larger, the discussed LIV test will 
most likely not be feasible with the current generation of IACTs. While the dependence of 
the critical LIV energy on the symmetry breaking scale is weak (see Sec. [TTJ) , for £i = 10 
E* reaches ~ 50 TeV. For a blazar such as those described in Fig. HJ it is apparent that 
the photons reaching Earth will be drastically diluted at energies below 50 TeV, so that the 
source will be completely unobservable at the energies where the effect of LIV would become 
substantial. Only if more nearby powerful extra-galactic sources are observed (see Fig. [2]), 
the gamma-rays at these high energies may be detected so that such weaker LIV scenarios 
can be examined. In any way, the examination of LIV scenarios with symmetry breaking 
around the Planck scale will be an important achievement, as such scenarios are commonly 
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predicted by theories and no generic test for them has yet been performed. 

V. DISCUSSION 

We show here that TeV blazar observations could distinguish LIV effects and provide 
unambiguous results. We can see in Fig. [2] that a source which displays an exponential 
cutoff due to EBL absorption in the presence of LIV at any energy below E* would show a 
re-emergence of gamma-rays at another energy above E*. The observation of such a strange 
phenomenon could clearly serve as a starting point leading us to LIV. Yet, because of the 
exotic nature of the LIV conjecture this would still not be convincing enough. However, if we 
observe the very-high-energy emission of a few different blazars, the ensemble of observations 
from different redshifts could serve to verify the LIV hypothesis. The different absorbed 
spectra must, of course, obey the same LIV parameters, hence we can assess all the re- 
emergence energies of different sources according to their cutoff energies. If a source cutoff 
is at an energy E% then we expect the re-emergence to be at the energy E2 that gives 
in dSJ) £thr{.E 2 ) = € t hr(Ex). Given an ensemble of cutoff and re-emergence energies, we can 
invoke a systematic statistical procedure, such as Maximum Likelihood, to estimate the most 
probable value of the LIV parameter, £. In case the above expectations are realized and the 
different sources display spectra that fit with high significance a certain LIV scenario, we 
will have very strong evidence for Lorentz violation, that would otherwise require extremely 
fine-tuned intrinsic effects or background spectrum. 

Alternatively, if we ever see an extra-galactic cutoff (or notice an increase in extra-galactic 
extinction) in a source's gamma-ray spectrum at an energy E, this will be a clear sign that 
a LIV model with E* < E is not a correct description of nature. We will have [see Eq. (j4])] 
a bound £ > [E n+2 (n + 1)(1 — 2~ n ) / {Am 2 e )] » /E p \. In the presence of LIV with symmetry 
breaking at the Planck scale, the optical depth never increases above E* = 23 TeV (note 
that far infra-red or CMB photons do not make any contribution to the optical depth). 
We must, of course, be convinced that the break in the spectrum is the result of extra- 
galactic attenuation and not an intrinsic feature of the source. This will be made possible by 
observations of cutoffs in the spectra of several different sources. Such observations will help 
us understand the intrinsic features of blazars, as breaks in the spectra that are not correlated 
with the sources' distances and appear at energies corresponding to the specifications of the 
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blazars could be interpreted as due to the blazars' mechanism, while breaks at energies 
that depend linearly on the sources' distances are most likely the result of cosmological 
absorption. A statistical analysis method can be used here as well to obtain a significant 
bound on the symmetry breaking scale. 

The notion of obtaining insight on LIV as described will of course remain a strictly 
theoretical idea, unless we can make the necessary observations and detect the radiation 
from the distant gamma-ray sources up to the relevant energies. Luckily, IACTs for very- 
high-energy gamma-ray astronomy are currently being rapidly enhanced. Two of the most 
important telescope arrays, MAGIC and HESS, are now being upgraded to Phase-II with 
the addition of improved larger telescopes. They will have an unprecedented light collection 
area. This will mainly extend the detection range to lower gamma-ray energies, but a 
slight increase in the detection rates of multi-TeV photons could, as discussed, be enough 
to make a significant discovery. A further improvement of the sensitivity in the multi-TeV 
gamma-ray band is expected with next-generation observatories such as CTA. Only nearby 
cosmological sources are relevant for the investigation of LIV (more distant sources would 
be attenuated at energies where LIV is negligible). Luckily again, at the source redshift of 
0.03, demonstrated in Figs. l3l4"l to be highly sensitive to LIV, there are already two blazars 
known to be powerful TeV gamma-ray emitters. By now Mkn-501 and Mkn-421 have been 
observed up to 20 TeV. With only a modest expectation from the upgraded IACTs, these or 
other nearby sources can be detected in the near future at slightly higher energies that will 
necessarily provide us input on the LIV possibility. The more distant gamma-ray sources 
will serve us for comparison with classical extra-galactic absorption features. 

Concluding, we have demonstrated a strategy that can test the LIV hypothesis by utiliz- 
ing its unique signature in observations of cosmological gamma-ray sources. The proposed 
examination of absorption features is independent of the background photon spectrum. We 
see that through a simple analysis of the observational data that is expected in the near 
future, we will either witness for the first time a clear sign of Lorentz violation or have a 
convincing generic argument against most theories that predict the break of the Lorentz 
symmetry at the Planck scale. 
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